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Abstract 

Recently we have demonstrated how to use partner symmetries for ob- 
taining noninvariant solutions of heavenly equations of Plebanski that 
govern heavenly gravitational metrics. In this paper, we present a class 
of scalar second-order PDEs with four variables, that possess partner 
symmetries and contain only second derivatives of the unknown. We 
present a general form of such a PDE together with recursion relations 
between partner symmetries. This general PDE is transformed to sev- 
eral simplest canonical forms containing the two heavenly equations 
of Plebanski among them and two other nonlinear equations which 
we call mixed heavenly equation and asymmetric heavenly equation. 
On an example of the mixed heavenly equation, we show how to use 
partner symmetries for obtaining noninvariant solutions of PDEs by 
a lift from invariant solutions. Finally, we present Ricci-flat self-dual 
metrics governed by solutions of the mixed heavenly equation and its 
Legendre transform. 



1 Introduction 



In his paper [1] , Plebanski introduced heavenly equations for a single poten- 
tial generating (anti-)self-dual heavenly metrics which satisfy complex vac- 
uum Einstein equations. Two real cross sections of these complex metrics, 
Kahler metrics with Euclidean or ultra-hyperbolic signature, are generated 
by the elliptic and hyperbolic complex Monge- Ampere equation [CM A) re- 
spectively, the particular cases of the first heavenly equation. Solutions of 
CM A play an important role in the theory of gravitational instantons [2], 
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where all gravitational metrics obtained so far, apart from the ones that we 
obtained lately [3,4], have Killing vectors, i.e. admit continuous symmetries. 
This implies symmetry reduction in the number of independent variables in 
metric components [5] , so these metrics actually live on manifolds of dimen- 
sions less than four. 

Recently we introduced the concept of partner symmetries and applied 
them to obtain noninvariant solutions of the complex Monge-Ampere equa- 
tion and the second heavenly equation of Plebanski [3,4,6,7]. Only such 
solutions could generate heavenly metrics with no Killing vectors, so that the 
metric components would depend on all four independent variables. Partner 
symmetries constitute a certain type of nonlocal symmetries and solutions 
invariant with respect to these nonlocal symmetries generically are solutions 
noninvariant in the usual sense, i.e. they depend on all four variables, so 
that no symmetry reduction in the number of independent variables occurs. 
The idea of using the invariance under nonlocal symmetries in order to get 
noninvariant solutions, suggested first by Dunajski and Mason [8,9], clarified 
for us the meaning of differential constraints, which we used earlier in [10,11] 
to derive non-invariant solutions of CM A. Since the partner symmetries 
and their use for lifting invariant solutions to noninvariant ones [4, 7] proved 
to be an appropriate tool for constructing noninvariant solutions of partial 
differential equations (PDEs) and a PDE taken at random would not admit 
partner symmetries, the natural question arises of how general this method 
of obtaining noninvariant solutions can be, or, in other words, what is the 
general form of equations that possess partner symmetries? 

To give a partial answer to this question, in this paper we present some 
results on a classification of the second-order PDEs of the general form 

-^(^tt) ^te) ^tyj ^tzj ^xxj '^xyj '^xzy ^yy^ ^yzj "^zz^ ^t; ^) ^) ■^^ U ^ ^) 

(1.1) 

that possess partner symmetries. Here u is the unknown that depends on 
the four independent variables t, x, y, z and the subscripts denote partial 
derivatives of m, e.g. Uu = d'^u/dt^, Utx = d'^u/dtdx . . .Although we have 
derived a complete set of equations for F such that equation (11.11) admits 
partner symmetries, we are currently able to give a general solution to these 
equations only for F that depends only on second derivatives of u. Thus, we 
obtain a classification of PDEs of the form 

Fi^Uft, Ufxj Ufy, Utzj Uxxi Uxyi Uxzi '^yy^ Uyz, Uzz) (1.2) 
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that possess partner symmetries. 

Our definition of partner symmetries requires tlie following two conditions 
to be satisfied: 

1. The symmetry condition for PDE (11. ip (determining equation for sym- 
metries) has the form of a two-dimensional divergence, that implies the 
existence of a single potential for each symmetry. 

2. The potential of each symmetry is itself a symmetry of (II. ip . i.e. a 
partner symmetry for the original symmetry. 

The first condition is satisfied in two steps: at first we require the symmetry 
condition to have the form of a four-dimensional divergence and then reduce 
this to a two-dimensional divergence by imposing additional constraints on F. 
We note that it were also possible to use a four-dimensional divergence form 
by introducing several potentials, as was shown, for example, by Bluman and 
Kumei [12], which would probably modify our concept of partner symmetries. 
This work is currently in progress. 

In section [2], we derive the restriction on the form of equation (11.11) implied 
by the requirement that the symmetry condition should have the form of a 
four- dimensional divergence: the left-hand side of the equation (11.11) itself 
should be a four-dimensional divergence, so that (II. ip becomes a conservation 
law. 

In section 121 we derive further conditions on F under which the four- 
dimensional divergence form of the symmetry condition is reduced to a two- 
dimensional divergence form which implies the existence of a single potential 
for each symmetry of (11.11) . 

In section HI we require that the potential of a symmetry should itself be 
a symmetry of the equation (11.11) and obtain the final set of equations for F. 
The definition of the symmetry potential then becomes a recursion relation 
for symmetries which generically maps any local symmetry into a certain 
nonlocal symmetry. To have an explicit form of this recursion relation, we still 
need a solution of the equations for F . We note that our symmetry potential 
is completely different from potential symmetries of Bluman and Kumei [12], 
where potentials are introduced not for symmetries but for PDEs, set in a 
divergence form, and symmetries are allowed to depend on these potential 
variables. 

In section [5], we attempt to solve the set of equations for F. The solution 
process in full generality turns out to be too lengthy and suggests many cases 
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and subcases to be considered. Therefore, here we restrict ourselves to the 
case where F in fll.ip depends only on the second derivatives of u and the 
equation takes the form (11.21) . Then we obtain a general solution for F in the 
left-hand side of (11.21) . that is, a general form of the second-order PDE with 
four variables containing only second derivatives of u that possesses partner 
symmetries, up to a change of notation for independent variables. We also 
obtain a recursion relation between partner symmetries in an explicit form. 

In section [6l we present a complete set of canonical forms, to which 
the general PDE with partner symmetries can be transformed by point and 
Legendre transformations, together with recursions for symmetries of these 
canonical equations. Among these canonical forms we find the first and sec- 
ond equations of Plebanski and two other nonlinear equations which we call 
mixed heavenly equation and asymmetric heavenly equation. The mixed 
heavenly equation turns out to be related by a partial Legendre transforma- 
tion to Husain's heavenly equation [13-15], which is an alternative form of the 
self-dual gravity equation related to the chiral model approach to self-dual 
gravity. 

In section [TJ we demonstrate an application of partner symmetries for 
finding noninvariant solutions of PDEs on an example of the mixed heavenly 
equation. We choose both symmetries in the recursion relations as transla- 
tional symmetries, with the recursions becoming differential constraints, and 
then show that Legendre transformation in two variables of both the equa- 
tion and two differential constraints leads to a set of three linear equations 
with constant coefficients. One of these equations depends only on three 
variables, containing the fourth variable merely as a parameter, and coin- 
cides with the Legendre transform of the translational symmetry reduction 
of the mixed heavenly equation, but expressed in new variables. Two other 
linear equations provide a lift of any solution of this equation, which is an 
invariant solution to the mixed heavenly equation, to a noninvariant solution 
that depends on all four variables. We present explicitly a linear combina- 
tion of exponential solutions and a polynomial solution as examples of such 
solutions. 

In section [HI we obtain Ricci-flat mixed heavenly metric in the self-dual 
gravity, governed by solutions of the mixed heavenly equation, by using a 
one-dimensional Legendre transformation of Husain's heavenly metric with 
a subsequent symmetrization of the transformed metric. Then we apply the 
linearizing Legendre transformation from section [7] to the mixed heavenly 
metric to obtain the Ricci-flat self-dual metric with a potential satisfying 



4 



the Legendre transformed mixed heavenly equation. We are now able to 
use our solutions of the latter equation, that are given in section [71 in the 
obtained metric or any other solutions of the above-mentioned three linear 
PDEs with constant coefficients. In this way we arrive at an explicit form 
of a Ricci-fiat self-dual metric with components generically depending on all 
four independent variables which, as a consequence, will admit no continuous 
symmetries. 

We have to mention also that a different classification of integrable three- 
and four-dimensional PDEs, that contain only second derivatives of the un- 
known, was given by Ferapontov et. al. in [16,17]. In this approach inte- 
grability is understood as the existence of sufficiently many hydrodynamic 
reductions, which is a requirement completely different from the existence of 
partner symmetries and therefore the results are also completely different. 

2 Divergence form of symmetry condition 

Let if he a. symmetry characteristic [18] of fll.ll) . Then the symmetry con- 
dition for symmetries ip admitted by (11. ip is determined by vanishing of the 
Frechet derivative of F on solutions of (11. ip 

A{ip) = Fuf + Fut^pt + Fu^Lf^ + Fuy^Py + Fu.ifz + Futt^u + 

+ FutyVty + Fut.^tz + Fu^:,Vxx + Fu^y^xy + Fu^^if^^ 

+ Fuyy^yy + Fuy^^yz + Fu^^^zz = 0, (2.1) 

where ipt = Dt^p, ipx = D^ip, . . . and Dt, D^, ■ ■ ■ denote operators of total 
derivatives with respect to t, x, . . . , e.g. 

A/ = df/dt + utdf/du + uttdf/dut + u^tdf/du^ + Uytdf/duy 
+ u^tdf/du^ + utttdf/dutt + Utxtdf/dutx + Utytdf/duty + ... 

After collecting all terms that can be written as total derivatives, the sym- 
metry condition (12. ip becomes 

A(M) + Dx{N) + Dy{L) + Dz{K) + Eu{F)^ = 0, (2.2) 

where Eu{F) denotes the Euler-Lagrange operator [18] applied to F 

E^{F) = D^iF.,J + DliF^J + DliF^J + Dl{F^J + D^F^J 
+ DtDy{F^,^) + DtDziF^J + DxDy{FuJ + AA(F„,J 
+ DyD,{F^J - A(i^«J - A.(i^.J - Dy{F^^) - D,{F^J + F„ (2.3) 
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and M, N,L,K are defined by 



+ 



F^^-D,{F^J--D.,{F^ 



1 



Dy{F^J--D,{F^J 



+ 



Fux ^x{Fy^^^) ^ Dfi^Fy^^^) ^ Dy(^Fy__^y) ^ Dz{Fy^^^ 



(2.4) 



+ 



111 

Fuy ~ F>y{^Fuyy) — — Dt{Fufy) " — Dx^yFu^y^ " " Dz^F^y^ 



+ 



if. 



The determining equation, transformed to the form (12.21) . takes the diver- 
gence form on solutions of ( 11. ID 

A(M) + D,(iV) + Dy{L) + Dz{K) = 0, (2.5) 

if and only if the Euler-Lagrange equation 

EuiF) = (2.6) 

is identically satisfied on solutions of F = 0, which is equivalent to the 
4-divergence form of the equation (11.11) itself [18]: 

F = Dt{P) + DxiQ) + Dy{R) + Dz{S) = 0, (2.7) 

where P,Q,R,S depend on the same set of variables as F in (II. ip . 



3 Two-dimensional divergence form of the 
symmetry condition 

In order to introduce a unique potential as a consequence of a symmetry 
condition, we have to convert the four- dimensional divergence in the left- 
hand side of the symmetry condition (12. 5p into a two-dimensional divergence. 
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say, in the variables t and x. To do this, we present L and K as the sum of 
total derivatives in t and x plus remaining terms which cannot be given in 
this form: 



1 



1 



1 



+ 



V, (3.1) 



K = Dtl -F^.^if ]+D,l -F^^^if + F^^^ip, + -F^y^^y 



+ 



F^^ - D^iFuJ - DtiFuJ - DxiF^J - -Dy{F^y 



(p. 



Using (13.11) in fl2.5p together with the definitions (12.41) and collecting terms 
with the total derivatives with respect to t and x, we convert (12.51) to the 
form 

+ [F^y - D,{F^J - DxiF^J] ^y + - D,{F^J - D^iF^J] 

+ {Dy [Fuy - Dy{F^J - D,{F^J - DxiF^J] (3.2) 

+ [F^^ - D^iF^J - D.iF^J - Dx{F^J] - DyD,{F^J]ip = 0, 

where M and N are defined by 



M = Fu,,(pt + ^Fu.^ipx + Fu.yipy + Fu,^(p, 



+ 
N 
+ 



F^, ~ D.iF^J ~ -DxiF^J 



Fu^^^x + ^Fu.^Lft + Fu,y(py + Fu,^(p, 



Fux F)x{Fu^^) 



(3.3) 



It is clear that in order to have the symmetry condition (13.21) to be a two- 
dimensional divergence in the variables t and x the coefficients of all the terms 
not included in the total derivatives Dt and should vanish on solutions 
of ([TT]): 

Fuyy = i^n,. = i^... = 0, (3.4) 
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F^^ - D,{F^J - D,{F^J = 0, - D,{F^J - D,{F^J = (3.5) 

whereas, as a consequence of fl3.4l) and (13. 5p . the coefficient of ip in (13. 2p 
vanishes identically and the symmetry condition (13. 2p becomes 

Dt{M)=D^{N). (3.6) 

Note that the symmetry condition and therefore all the equations (13. 2p . (13. 4p . 
(13. Sp . and (13. 6p should be satisfied not identically but only on solutions of 
the original PDE (II. ip and hence they should be (differential) consequences 
of F = 0. 

Condition (13.60 is equivalent to the local existence of the potential 
defined by 

^, = N = N + At, ^, = M=M + A,, A = ujip, (3.7) 

where u may depend on t, x, y, z, u and the first and second derivatives of u. 
Here the terms with the derivatives of A are added in order to have the most 
general definition of the potential if). Now, the symmetry condition (13. 2p can 
be written as 

Dt{M)=D^{N) (3.8) 

on solutions of F = 0. 

4 Existence conditions for partner 
symmetries 

Our second requirement is that the potential ip should also be a symmetry 
of the PDE (II. ip . i.e. a partner symmetry for the original symmetry ip, so 
that (13. 7p becomes a recursion relation for symmetries. Then the symmetry 
condition in the two-dimensional divergence form (13. 8p with ip replaced by 
Ip, defined by (13.70 . should be satisfied on solutions of the equation (11.10 

a(m) =/^. (iv)+F, (4.1) 
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where M and N are obtained from M and respectively by replacing 
with ip in fl3.3l) 



+ 
iV 

+ 



Fut- F>t{Fuu) - 2^xiFutJ 



(4.2) 



Fu^ - D^{Fu^J - -Dt{Fut^ 



The term F has the form 

F = ^iDtiF) + uD^iF) + pDy{F) + XD,{F) + aF 



(4.3) 



and it accounts for the fact that equation (14.11) should be satisfied only on 
solutions of (11.11) (a consequence of proposition 2.10 in [18], similar to formula 
(2.26) therein). Terms with ipy, ip^, and in (14.11) cannot be balanced by 
any other terms and therefore they should vanish separately on solutions of 
f 1 1.11) yielding 



and 



D,{F^J + D,{F.) = py, A(Fn.J + D,{F^^ 



D,{F^,) + D.(F„J - Dl{F^,,) - DliF^J - D^F^ 



(4.4) 



(4.5) 



respectively, where the terms F^ and F^ are of the same form (14. 3 p but 
with different coefficients /i, z/, p, A, a. Equations (14.41) together with (13.51) 
and equation (14.51) together with (12.61) imply 



Fu^ = 0, =0, F„ = 0. 



(4.6) 



In all other terms in (14.11) , we replace ipt and ipx by the expressions (13. 7p . We 
note that, due to the definition (13. 7p of the potential ip and its consequence 
(13. 8p (equivalent to (13.60 ). ip satisfies the symmetry condition (12. ip . which 
cancels all the terms proportional to u in (14. ip . All other terms in (14. ip 
with second derivatives of cp are cancelled identically. The remaining terms 
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are proportional to Lpt, ^x-, <^y, <^z, and (/?, so that these five groups of terms 
should vanish separately on solutions of the equation F = to give the 
following five equations respectively: 

- ^ [FuMF^J + F^,^D,{F^J] + 2F^M^) + FuM^j) 

+ F^,^Dy{uj) + F^^M^) = A, (4.7) 



-\ F>t{FuttFuxx) - 2 F'tiF^tJ + FutyDy{Fu^^) + Fut,D^{Fu^J 



^ [F^^^DyiF^J + F^^D,{F^J] - 2F^^^D,{uj) - F^^^D,{u) 



l} = A, (4.8) 



Dxi^FutyFuxx) 2 -^^^Futx-Fuxy) t^Fy^^^ F^^^) + ^ Dti^Fuf^Fufy 
"I" FuxyDyi^Fy-ty) + Fu^^Dz{Futy) — Fui-yDyiyFuxy) ~ Fui-^D zi^Fy^^y) 
+ F^,yDt{uj) + FuxyD^iuj) = Fa, (4.9) 

Dx{,FutzFuxx) ~ 2 -^^i-Futx-Fuxz) ~ ■F)t{Fu^^Fu^^^ + — F)f[Fu^^Fu^^) 

~l~ FuxyDyiFu^^) + Fu^^D z(yFu^^) — F^yDyi^Fu^^) — Fu^^Dz{Fu^^) 

+ A(a;) + Fy^M^) = (4.10) 



F„,B - F„^A + F„„ A(5) - i^«..^x(^) + ^ Fy^^ [D,{B) - A(^)] 

+F„,^F>,(S) + F^^D,{B) - Fy^yDyiA) - F^^D^) + i(u;) = Fq, 

(4.11) 

where 

A = A(i^«J + ^^.(i^„.J-i^«„ s = A(i^„,J + ^A(i^«.J-i^«. (4.12) 
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and A is the operator of the symmetry condition (I2.1I] . Here the terms F,, 
of the form (14 ■3p but with different coefficients, account for the fact that the 
equations should be satisfied only on solutions of (11. ip . 

We note that in the notation (14.120 equation (14. 5p simplifies to 

Dt{A) + D^{B) = F. (4.13) 

5 Equations that admit partner symmetries 
and recursion relation for symmetries 

We proceed now to solve the existence conditions for partner symmetries 
dSD, dSD, (USD, dUD, dUni), (BHD and dUSD for the unknown left-hand 
side F of the equation (11.11) and u in the definition (13. 7p of the potential ip. 
We split these equations in third derivatives of u obtaining over-determined 
sets of equations which can be easily solved. Our strategy is to choose the 
function u and the coefficients fi,^, p,X,a in the terms F of the form (14.31) in 
such a way as to have minimum restrictions on the form F of equation (II. ip . 

We start with the equations (14.40 since they do not contain u. Our 
strategy results in vanishing of /i, u, p, A and a in F^ and F^ that implies the 
linear dependence of F on Uty,Uxy,Utz and Uxz, so that the solution of the 
equations (14. 4p has the form 

F = ai{y, z){utyUxz - UtzU^y) + a2{utxUty - UuUxy) + as{utyUxx - UtxUxy) 

+ a4{UtxUtz - UttUxz) + a^iUtzUxx - UtxUxz) + hUxy + b2Uty 

+ huxz + biUtz + g:i{utu Utx, Uxx, Ut, Ux, t, X, y, z), (5.1) 

where the coefficients 02, 03, 04, a^, bi, 62, &3 and 64 are functions of Ut, Ux, t, x, 
y, z that satisfy a certain overdetermined set of partial differential equations. 

To simplify the analysis, we assume from now on that all the coefficients in 
(15. ip are constants, so that all the equations for the coefficients are identically 
satisfied, and depends only on the second derivatives Utt,Utx,Uxx- As a 
consequence, the left-hand side F of our equation (11.10 depends only on 
second derivatives of u and it takes the form (II. 2p . 

With these restrictions, we substitute the expression (15.10 for F in the 
remaining six equations (gH), (glD, (glD, KTOh . flCTjl and KT3^ . The 
resulting equations are split in third derivatives of u into over-determined sets 
of equations, where we choose the function u and the coefficients fi,^, p, X,a 
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in the terms of the form fl4.3l) in such a way as to obtain minimum restrictions 
on the form of F . It turns out that all these six equations determine only 
the form of the function g^^ 

gsiutt, utx, u^x) = deiuttUxx - Utx) + hutt + '^hutx + hu.^^ + 60, (5.2) 
so that the equation (11.11) becomes 

F = aiiUtyU^z - UtzU^y) + a2{UtxUty - UttU^y) + a^iUtyUxx - UtxUxy) 

+ a^iuixUtz - UttUxz) + a^iutzU^x - UtxUxz) + aeiuttUxx - ^L) 

+ hu^y + b2Uty + hu^z + hutz + hutt + 2h^Uta: + hu^^ + 60 = (5.3) 

together with the following solution for uj 

= {a2Uty + tt'iUxy + a^utz + a^u^z) + ^0, (5.4) 

where all the coefficients are constants. Using (15. 4p in the equations (13. 7p . 
that define the symmetry potential in terms of the symmetry y9, we obtain 
the recursion relation between partner symmetries of the equation (15.31) 

A = -{a2Uty + a^z - deUtx + - uJo)y:>t 

- [a^uty + a^utz + a^utt + b-j)ipx + [aiUtz + a2Utt + a^Utx - bi)(py 

+ ( - aiUty + a^Utt + a^Utx - b^)'^z, (5.5) 
"^x = -{a'2Uxy + a^Uxz - cigUxx - ^sj^^t 

- {asUxy + a-bUxz + anutx - be- loq) (fx 

+ {aiUxz + «2Mte + cisUxx + b2)^y + { - aiU^y + a^utx + a^Uxx + b/^^z- 

Here, by construction, both ip and ip satisfy the symmetry condition (12.11) in 
the divergence form (13.81) and (14. ip respectively, on solutions of (11.21) . and 
hence the transformation (15.51) maps any symmetry ip of the equation (15.31) 
again into its symmetry ip. 

6 Canonical forms of the PDEs that admit 
partner symmetries 

Due to the random choice of original variables, both the form of the equa- 
tion (15. 3p . that admits partner symmetries, and the recursion relation (15.51) 
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contain false generality. Therefore, we present here simple canonical forms, 
to which the equation fl5.3p can be transformed by point and Legendre trans- 
formations, and the corresponding recursion relations for symmetries. This 
will also make up for our casual choice of variables t and x for the two- 
dimensional divergence form. The transformations, providing the proof of 
the results presented here, will be published elsewhere. 

Case I: ai 7^ 0. 

In this case we can make ai = 1 by dividing (15. Sp over ai. The equation (15. 3p 
can be reduced to the form 

F = UtyU^z - UtzUxy + ^{UttUxx " ufx) + + Bu^^ + CUtx 

+ hiu^y + h2Uty + huxz + hutz + 60 = 0. (6.1) 
Consider here the following subcases. 
Subcase la: F = 0. 

Then, by the change of the unknown the equation (16. ip becomes 

UtyU^z - UtzUj^y + Autt + Bu.^^ + Cutx + -D = 0. (6.2) 
Subcase lal: A = B = C = D = 

Then equation (16. 2p reduces to the homogeneous version of the first heavenly 
equation of Plebanski 

UtyUxz - UtzUxy = 0. (6.3) 

Subcase Ia2: A = B = C = 0, D^O. 

Then we can set D = — 1 and the equation (16.20 becomes the first heavenly 
equation of Plebanski 

UtyUxz - UtzUxy = 1. (6.4) 

In the cases lal and Ia2, the recursion relation for symmetries (15. 5 p becomes 

Ipt = UJQl-Pt + Utz<-Py - UtyiPz, Ipx = i^O<^x +Uxz<-Py - UxyiPz- (6.5) 

Subcase Ia3: {A, B,C) ^ (0,0,0). 

Then we can always make D 0. By a combination of Legendre and point 
transformations the equation (16.20 takes the canonical form 

UtyUxz - UtzUxy + UttUxx - = (6-6) 
where e = ±1. We call (16. 6p the mixed heavenly equation. 
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The homogeneous version of the mixed heavenly equation fl6.6p . with 
e = 0, can be transformed to the first heavenly equation by an appropriate 
Legendre transformation. 

The recursion relation fl5.5p for symmetries of equation fl6.6p becomes 

A = {Utx + ^o)Vt - Utt^x + Utz^Py - Utyifiz, 

■ipx = UxxVt - {Utx - i^o)<-Px + Uxz<-Py " U^y<^z- (6.7) 

Recently we became aware of the relation of the mixed heavenly equation 
to the Husain's heavenly equation (at e = +1) [13, 14] arising in the chiral 
model approach to self-dual gravity 

VtyVpz - VtzVpy + vtt + evpp = 0. (6.8) 

Husain's equation can be obtained from the mixed heavenly equation by the 
partial Legendre transformation in x 

p = Ux, v{t,p,y, z) = u — xUx, X = —Vp, u = V — pVp. (6.9) 

We note that (16. 8p could also be obtained as a canonical equation in the 
subcase Ia3 of the general equation (15.31) with the replacement u ^ v,x ^ p 
with the following choice of the coefficients in (15. 3p : Oi = 1, 65 = 1, 67 = £ 
and all other coefficients vanishing. Then from (15. 5p . with this change of 
notation, we obtain the recursion for partner symmetries of equation (16. Sp 

i/jt = uJo(pt+VtziPy-Vtyipz-£^p, = ^t+u:o(pp + Vpz(py-Vpy(pz- (6.10) 

Subcase lb: T ^ 0. 

In this case, the equation (16.11) can be transformed to the same equation 
O- 

Case II: ai = 0. 

Subcase Ila: 02 = as = ^4 = 05 = 0, ag 7^ 0. 

Then we can set = 1 and the equation (15. 3p becomes 

F = UttUxx - ul^ + b^utt + 266Mte + huxx 
+ hiu^y + h2Uty + hu^z + hutz + 60 = 0. (6.11) 
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The equation 06.111) can be transformed to the equation 

UttUxx - + blUxy + b2Uty + huxz + hutz = 0. (6.12) 

We consider only the case when (61,62) 7^ (0,0) and (63,64) 7^ (0,0). Other- 
wise the equation (16.121) will determine a function of less than four variables. 

Case Hal: 6164 - 6263 ^ 0. 

Then by a change of variables the equation (I6.12p takes the form of the second 
heavenly equation of Plebanski 

UttUxx - Utx + "^xy + Utz =0- (6.13) 

The recursion flS.Sp for symmetries of (16.131) takes the form 

Ipt = {Utx + UJo)(pt - Utt^x - ^y, Ipx = Ua,x^t - {Utx - i^o)Vx + Vz, (6.14) 

which at ciJo = coincides, up to the change ip t— > —ip, with our previous 
result [6]. 

Case IIa2: 6164 - 6263 = 0. 

Then, by choosing a certain linear combination of y and z, we obtain the 
equation which determines a function of only three variables. 

Subcase lib: ai = a2 = = = = = 0. 
Then (15. 3p reduces to the hnear equation 

b5Utt + 2bQUtx + huxx + hiUn^y + h2Uty + hu.j,z + h^utz + 60 = 0. (6.15) 

The recursion (15. 5p for symmetries of (16.150 becomes 

= -(65 - t^o)v't - hVx - bi(py - h(pz, 

i^x = h^t + (65 + ^o)^x + 62<^2y + h^ifz- (6.16) 

Subcase He: ai = 0, (02,03,04,05) 7^ (0,0,0,0). 

The equation (15. 3p with oi = can be reduced to the canonical form 

UtxUty - UttUxy + autz + huxz + cuxx = 0, (6.17) 

up to a possible change of notation for independent variables x ^ t and/or 
z y. Here we can set o = 1 by an appropriate scaling of variables. We 
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call this equation asymmetric heavenly equation. At 6 = it becomes the 
so-called evolution form of the second heavenly equation [15,19,20]. 
The recursion relation fl5.5p for symmetries of (16.171) takes the form 

-ipt = -{Uty - i^o)Vt - Cipx + Utt^y - b^pz 

i>x = -Uxy^t + ^^O^x + Utx^y + a^z- (6.18) 

7 Lift from invariant to noninvariant 

solutions of the mixed heavenly equation 

In the method of partner symmetries, we consider a nonlocal symmetry with 
the characteristic fj = (p — Rep, where (p is any point symmetry of our equa- 
tion and R is the recursion operator determined by recursion relations (15. 5p . 
generating a nonlocal symmetry ijj = Rip from a point symmetry ip. We 
search for solutions invariant with respect to a nonlocal symmetry fj, deter- 
mined by the condition ip — Rip = [3], so that we can obtain this invariance 
condition by formally replacing ip by a point symmetry ip : ip = ip in the 
recursion relations (15.51) . This does not mean symmetry reduction, so that 
generically these solutions depend on all four variables and so they are still 
noninvariant solutions in the usual sense. Contact symmetries can also be 
used for ip and/or ip. 

Here we demonstrate the application of partner symmetries for obtaining 
noninvariant solutions of canonical PDEs and, in particular, a lift from in- 
variant to noninvariant solutions. We choose mixed heavenly equation (16.61) 
as an example, possessing the recursion for symmetries (16. 7p . where we set 
i^Q = 0. The equation (16.61) admits the obvious translational symmetry with 
the generator X = dx + dz. 

Solutions, invariant under this symmetry, have the form u = u{s,t,y), 
where s = x — z, since they do not change under the simultaneous shift in x 
and z. Then u satisfies the reduced equation 

UtsUsy — UtyUss + UftUss ~ Ufs — ^i 

obtained from (16.60 by the symmetry reduction. Under the Legendre trans- 
formation 

r = Us, v{r,t,y) = u — sUs, s = —Vr, u = v — rvr (7.2) 
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the equation fl7.ip is linearized in the form 



Vtt + eVrr - Vty = 0. (7.3) 

By using partner symmetries, we shall show that solutions of the linear equa- 
tion (17.31) ■ i.e. invariant solutions of Legendre transformed mixed heavenly 
equation, being written in certain new coordinates, can be lifted up to non- 
invariant solutions of the latter equation. 

As was explained at the beginning of this section, we formally replace ip, 
that is generated from a point symmetry in the recursion relations (16.71) 
(with cuq = 0), by a point symmetry ip. Here we choose both ip and ip = (p 
to be the indicated above combination of translations in x and z with the 
characteristic ip = tp = Ux + Uz, so that (16. 7p becomes 

Uxx ~\~ Vjxz UtzUxx UtxUxz ~l~ UxzUyz UxyUzzy C^'^) 
Utx + Utz = Utx{Utx + Utz) — Utt{Uxx + Uxz) 

+ Utz{Uxy + Uyz) - UtyiUxz + U^z)- (7.5) 

With the aid of (16.61) . the equation (17.51) takes the form 

Utx + Utz = UtxUtz - UttUxz + UtzUyz - UtyU^z " £■ (7.6) 

After the Legendre transformation 
p = Ux, q = Uz, v{p,q,t,y) = u-xux-zuz, x = -Vp, z = -Vg{7.7) 
the equations (17. 4p and (17.61) take the form 

Vpq = Vqq + Vtq - Vpy, (7.8) 
Vpq{Vtq - eVpq + Vtt) = Vpp{-eVqq + Vtq + Vty), (7.9) 

where the equation (17. 8p was used in the Legendre transform of (17. 6p to arrive 
at (17. 9p . The equation (17. 9p can be set into a linear form 

XVpq = Vtq - eVqg + Vty, (7.10) 
XVpp = Vtq - eVpq + Vtt (7.11) 

by introducing an extra unknown A depending on all the variables. Solving 
algebraically the system of the three linear equations (17. Sp . (I7.10p and (17.110 



17 



with respect to the principal derivatives Vty^ Vtq and Vpy in terms of the 
remaining parametric derivatives in the form 

Vty = £{Vqq " Vpq) + \{Vpq - Vpp) + Vft, (7.12) 

Vtq = evpq + Xvpp - vtt, (7-13) 

Vpy = (e - l)Vpq + XVpp + Vqq - Vtt, (7.14) 

we easily check that all cross derivatives of the left-hand sides coincide as a 
consequence of these equations, so that this system of PDEs does not have 
nontrivial integrability conditions. The mixed heavenly equation (16.61) after 
the Legendre transformation (17. 7p becomes 

VtqVpy - VpqVty + VttVqq - V^^ + e{VppVqq - V^^) = 0. (7.15) 

The Legendre transformed mixed heavenly equation (17.151) obviously consti- 
tutes another particular case of our general equation (15. 3p . up to a change 
of notation of the dependent and independent variables. The linear equa- 
tions (I7A2D . (I7:T3|) . dHH) together with (rrT5|) imply that A = -e as far as 
VppVqq — Vpy 7^ 0. Then the Legendre transformed mixed heavenly equation 
(I7.15P becomes an algebraic consequence of these three linear PDEs with 
constant coefficients. 

In the case when e = —1 and hence A = 1, the equations (17. 8p and (I7.1UI) 
imply 

Vty + Vpy = 0, (7.16) 
which can be integrated to yield the linear first-order equation 

Vt + Vp = C{t,p,q). (7.17) 

This obviously leads to dependence of v on the characteristic combination 
t—p and thus determines invariant solutions. In this case we have a symmetry 
reduction and no lift to noninvariant solutions. 

In the case when e = 1 and hence A = — 1, the equations (17.80 . (17.100 
and (17. lip do not imply any linear first-order consequences, so there is no 
symmetry reduction of the number of variables in this case and invariant 
solutions generically do not arise. Under the change of variables {q,p,t) i-^ 
{q,ri = p + t,^ = p — t), the equation (I7.13P takes the form of the linear 
reduced equation (17.30 but written in the new variables r], ^ and q 

Vrjn + V^^ - V^q = (7.18) 



18 



and containing the fourth argument y of the unknown v as a. parameter. Any 
solution of this hnear equation depends on the three variables r;, ^ and q with 
the "constants" of integration depending on the fourth variable y. Certain 
appropriate linearly independent combinations of two other equations (17.121) 
and (I7.14p . with the use of (17.181) . in the new variables take the form 

V^q- V^q+ V^y = (7.19) 
V£_q + V^q - Vqq + V^y = 0. (7.20) 

These two equations determine the ^/-dependence of the "constants" of in- 
tegration in the solution of (I7.18P and hence we obtain the lift of invariant 
solutions of the Legendre transformed mixed heavenly equation (17.151) to 
noninvariant solutions of this equation. 

It is easy to obtain an infinite set of exact solutions to linear equations 
with constant coefficients. Indeed, we can try the exponential dependence of 
V on 7] 

V = exp {ar] + b{y))f{^,q,y), 

so that (17.181) becomes 

/«?-/«, + aV = 0. (7.21) 
For the solution of this equation we can try the following ansatz 

/ = A cos («e + Pq + 0{y)) + B sin («e + Pq + e{y)). (7.22) 



The expression (17.221) satisfies (I7.2ip only if a = -La^a — (3, so that 
V = exp ^ ± a^/a — (3 r] + b{y)^ 

X A cos {a^ + pq + 9{y)) + B sin {a^ + Pq + 9{y)) . (7.23) 
The expression (17.230 satisfies (I7.19P if 



e{y) = -(3y, b{y) = ±f3d^y 



a 



so that (17.230 finally becomes 



exp ( ±^/a{a - (3) {r] + -y 



a 



X I A cos K + Piq -y)] + B sin K + f^iq -?/)]}• (7-24) 
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Surprisingly enough, this expression satisfies identically the third equation 
fl7.20p . though it is linearly independent of the other two equations. 

Any linear combination of the solutions of the form fl7.24p is again a 
solution of the three linear equations (17.181) fl7.19p and (17.200 and hence, 
with rj = p + t, = p — t, it will satisfy the nonlinear Legendre transformed 
mixed heavenly equation (I7.15P at e = 1, since it is a consequence of these 
linear equations. In the case of a discrete spectrum, we can choose, for 
example, the following linear combination 



^exp f ±v^ 



ai{ai - /3i) [7] + —y 



Oti 



+ 5iSin [aii + (3i{q - y)] 



Ai cos [aii + (3i{q - y) 

(7.25) 



where a,, Ai and Bi stand for arbitrary constants. This is an example 
of a solution to (ETSj), which is obviously noninvariant because it clearly 
depends on four independent combinations of the variables r], ^, q and y. For 
the case of a continuous spectrum, the sum in (17.251) should be replaced by 
an integral. 

There is also a class of polynomial solutions. We start with the ansatz 



V = A{ri, y)y + 5(r/, ^, y)q + C{ri, ^, y). 



(7.26) 



The expression (I7.26P will satisfy linear equations (17.181) (17.191) and (I7.20p if 
the coefficients have the form 



A{7],i,y) = ?, 
B{7],^,y) = 3 



Hv'-f) + '2hr]^ + ky' 



{4g + h){e-v') + 2m-h)r]^ 



C{r], y) = k7]y' + 3 \h{r]' - e) - Sgv^] + fi^V 



(7.27) 



he + 8gv' + I2gr,^i - 3(4(7 + h)r,e + 3/i(r/2 - e) - Qf^V^ 



y, 



where /, g, h, k and /i are arbitrary constants. 

For solutions independent of 77 we obtain, for example, v = + q — y)^ 
and, since all the three equations are linear, the sum of this solution and 
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fl7.27p is again a solution, so that 



V = A{i^, i,y)\ + B{7^, y)q + Civ, y) + D{C + q- y)\ (7.28) 

with A, B and C defined by f l7.27p and constant D, will satisfy the non- 
linear Legendre transformed mixed heavenly equation fl7.15p . This solution 
generically depends on all the four variables and hence it is a noninvariant 
solution, i.e. it does not admit Lie symmetries. More general polynomial 
solutions can easily be constructed. The sum of the exponential solution 
f l7.25p and a polynomial solution again satisfies (17.150 . 

8 Ricci-flat metrics governed by the mixed 
heavenly equation and its Legendre 
transform 

In section [TJ we have obtained noninvariant solutions (I7.25P and (I7.28P of 
the Legendre transformed mixed heavenly equation ( 17.150 . In order to get 
the corresponding solution of the mixed heavenly equation (16. 6p . we had to 
perform the Legendre transformation of solutions (17.250 and (17.280 . inverse 
to (17.70 . which is quite a difficult problem. 

Instead, we shall proceed, as we did before in [3,4,6], by taking into 
account that, similar to the complex Monge- Ampere equation and second 
heavenly equation of Plebanski, the mixed heavenly equation determines a 
potential that governs Ricci-fiat metrics in the self-dual gravity. If we are 
interested only in such metrics as our final result, then instead of perform- 
ing the inverse Legendre transformation of our solution, we make the direct 
Legendre transformation (17. 7p of the metric related to the mixed heavenly 
equation. Then our solutions (17.250 and (17.280 of the Legendre transformed 
mixed heavenly equation (I7.15P at e = 1, or any other its solutions deter- 
mined by the linear equations (17.190 . (17.180 and (I7.20p . will yield a potential 
governing the Legendre transformed mixed heavenly metric. 

In order to obtain Ricci-fiat metrics related to the mixed heavenly equa- 
tion, we start with Husain's heavenly metric and then use the relation be- 
tween the Husain's equation and mixed heavenly equation. Husain's heavenly 
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metric has the form [13] 



ds = 2 ( ujtdt + Updp + 

where 



ujt = ^tydy + Atzdz, Up = Apydy + Ap^dz, Atp = AtyAp^ - AtzApy 

with A{t,p, y, z) satisfying the Husain's equation 

Att + App + AtzApy - AtyApz = 0. (8.2) 

By the one-dimensional Legendre transformation 

A = u — xux, p = —Ux, X = Ap, u = A — pAp, (8.3) 

where the inverse transformation is also given, the Husain's equation (18. 2p is 
mapped into the mixed heavenly equation with e = +1 

UtyUxz - UtzUxy + UttUxx - = 1 (8-4) 

for the unknown u{t,x,y, z). Performing the transformation (18. 3p of the 
Husain's metric (18.11) . we obtain the metric governed by equation (18.41) 

ds'^ = 2 lujtdt + ujxdx H ^ \{uxxi^t - Utx^^xf + (A + l)ul] \ , (8.5) 

L UxxA J 

where 

Ut = Utydy + Utzdz, Ux = Uxydy + Uxzdz, A = UtzUxy - UtyUxz- (8.6) 

By using a REDUCE program, it has been checked that the metric (18. 5p is 
Ricci-fiat as a consequence of equation (18.41) . 

The asymmetry of the metric (18.50 in variables t and x is caused by 
the Legendre transformation (18.30 between p and x, which leaves t untrans- 
formed. To amend this lack of symmetry, we symmetrize the metric (18. 5p in 
t ^ X and y z and then introduce t ± x and y ± z as new coordinates, 
which we call again t, x, y, z. The resulting mixed heavenly metric has the 
form 

ds^ = 2 i^utdt + Uxdx + ^ {uxx^^t " ^Utx^JtUJx + UuujI) | . (8.7) 
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This metric is also Ricci-flat, provided that the potential u satisfies the mixed 
heavenly equation (16.61) for both signs of e. Now we apply the Legendre 
transformation (17. 7p to the mixed heavenly metric (18.71) with the result 



where 5 = VtyVpq — VtqVpy and the metric potential v{t,p,q,y) should satisfy 
the Legendre transformed mixed heavenly equation (17.151) . With the latter 
condition satisfied, by using REDUCE we have checked that the metric (18.81) 
is Ricci-fiat and calculated the Riemann curvature tensor components for 
an arbitrary v satisfying (17.151) . The expressions for these components are 
too lengthy to be presented for publication. However, the denominators of 
the Riemann tensor components are simple, so that possible singularities 
of the curvature tensor either coincide with the singularities of the metric 
(18.81) . being at S = VtyVpq — VtqVpy = 0, or are located at Vqq = 0, for v 
being a linear function of q. For the polynomial solution (17.281) the condition 
6 = could be satisfied only if all the essential coefficients in (I7.28P vanished: 
h = g = fi = D = 0, which would contradict the non-invariance of this 
solution. The only singularity of the metric corresponding to (I7.28P is located 
at infinity. 

As it was shown in section [71 we can use any solution of the three lin- 
ear equations (17.180 (I7.19P and (I7.20p . which imply (I7.15P at e = 1 as their 
algebraic consequence. In particular, we can use the noninvariant solutions 
(I7.25P and (I7.28P for v in the metric (18. 8p . For noninvariant solutions, there 
will be no symmetry reduction, so that v will depend on all the four indepen- 
dent variables, which is a necessary (and often sufficient) condition for the 
metric components in (18. 8p to depend also on all the four independent vari- 
ables. For the exponential solutions to the complex Monge- Ampere equation 
and to the second heavenly equation, which are similar to (I7.25p . we have 
proved in [6] that the corresponding Kahler metric and the second heavenly 
metric admit no Killing vectors. Similarly, for the solution (17.250 we also 
expect that the Legendre transformed mixed heavenly metric (18.80 will ad- 
mit no Killing vectors and hence no symmetry reduction in the number of 
independent variables will occur. 



ds^ = 2 | - {vtt + evpp){vqtdt + Vqpdp + Vqqdqf 
\ 2 

+ {Vqtdt + Vqpdp + Vqqdq) -dq + - {VtqVty + eVpqVpy)dy 

+ Vytdt + Vypdp + -^{Vyt + eVyp)dy dy> 
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9 Conclusion 



In the theory of gravitational instantons, heavenly metrics with no Killing 
vectors (no continuous symmetries) can only be generated by noninvariant 
solutions of CM A. Therefore, we are faced with the problem of obtaining 
noninvariant solutions of partial differential equations. Partner symmetries 
proved to be an appropriate tool for solving such a problem because non- 
invariant solutions can be obtained as solutions invariant with respect to 
a certain nonlocal symmetry closely related to partner symmetries. Thus, 
the existence of partner symmetries for a given PDE is necessary to apply 
this method. In this paper, we have obtained a general form of the scalar 
second-order PDE in four variables, containing only second derivatives of 
the unknown, that possesses partner symmetries. Using point and Legendre 
transformations, we have transformed this general equation to different sim- 
plest canonical forms and so presented a classification of inequivalent equa- 
tions which admit partner symmetries, together with recursion relations for 
symmetries. Among these equations we find the well-known first and second 
heavenly equations of Plebafiski and two other nonlinear equations which 
we have called mixed heavenly equation and asymmetric heavenly equation. 
The mixed heavenly equation is related by a partial Legendre transforma- 
tion to Husain's heavenly equation arising in the chiral model approach to 
self-dual gravity. A particular case of the asymmetric heavenly equation is 
the evolution form of the second heavenly equation. 

We ignored here all the cases when the canonical equation explicitly con- 
tains only three variables. We leave for the future a classification of PDEs 
with three variables, that admit partner symmetries. 

As an example of application of partner symmetries, we have shown how 
to construct noninvariant solutions of the Legendre transformed mixed heav- 
enly equation. By applying Legendre transformation in two variables, the 
latter equation and differential constraints, that are obtained from recur- 
sion relations for partner symmetries, have been transformed to a set of 
three linear equations with constant coefficients, that imply the Legendre 
transformed mixed heavenly equation as their algebraic consequence. One of 
these equations involves only three variables and formally coincides with a 
certain reduced equation, which determines invariant solutions of the Legen- 
dre transformed mixed heavenly equation, but written in new variables and 
containing also the fourth variable as a parameter. Two other equations, 
involving all the four variables, provide a lift from invariant to noninvariant 
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solutions of the Legendre transformed mixed heavenly equation. We have 
obtained Ricci-flat metrics governed by the mixed heavenly equation and the 
Legendre transformed mixed heavenly equation. Using any noninvariant so- 
lution of the three linear PDEs, we satisfy the necessary condition of arriving 
at Ricci-flat metrics with metric components depending on all four indepen- 
dent variables. Such metrics will admit no continuous symmetries and no 
Killing vectors. 

Thus, we conclude that, for a scalar second-order PDE with four inde- 
pendent variables, the existence of partner symmetries happens to be a char- 
acteristic feature of the equations that describe self-dual gravity in different 
variables. The partner symmetries provide a tool for obtaining noninvariant 
solutions of these equations and Ricci-flat self-dual metrics with no Killing 
vectors. 
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